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phase transitions. Results obtained from numerical simulations corroborate the predictions of the 
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I. INTRODUCTION 

The microcanonical analysis of phase transitions has 
gained growing interest in recent years due to the in- 
equivalence of the canonical and the microcanonical de- 
scriptions of finite systems with short-range interac- 
tions as well as of systems with lon g-ra nge interactions 
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Whereas the discrimination of the order of a phase tran- 
sition from finite-size data can be a difficult task in 
the canonical ensemble, microcanonical quantities usu- 
ally yield clear signatures which reveal the continuous or 
the discontinuous character of a phase transition 0. In 
small microcanonical systems which exhibit a continuous 
transition in the thermodynamic limit the character of 
the transition is signaled by typical features of symme- 
try breaking, as e. g. the abrupt onset of a non-zero or- 
der parameter when a pseudo-critical point is approached 
from above. These features are accompanied by sing ular 
physical quantities [HlllliHilllllillilEE El 113 . 
The associated singularities, however, are not the con- 
sequence of a non-analytic microcanonical entropy and 
can therefore be characterized by classical critical expo- 
nents [Hill. Similarly, intri guing features are also re- 
vealed in the microcanonical entropy of small systems 
with a discontinuous transition in the infinite volume 
limit aHHHHSIHIlllll. A typical back-bending of 
the microcanonical caloric curve is observed, leading to 
a negative heat capacity. This property which has also 
been observed in recent experiments |3 113 j has been 
used as a tool for the determination of the first order 
character of phase transitions in microcanonical systems. 
In the canonical context the first order character is then 
signaled by a double-peak structure in certain histograms 
(e.g. mlsllsi)- 

However, some systems (as for example the Baxter- 
Wu model ) , with a continuous phase transition 
in the thermodynamic limit, have been reported to show 
a similar double-peak structure in canonical histograms 



or, equivalently, a negative microcanonical specific heat 
|4lL 14^ . It follows from this observation that the mere 
existence of a convex dip in the entropy is not sufficient 
for the identification of a discontinuous transition ■ 

In this work we examine more closely the finite-size 
behaviour of the microcanonical specific heat and of the 
convex dip in the entropy function of this kind of systems, 
thereby taking the Baxter- Wu model and the four-state 
Potts model in two dimensions as examples. In partic- 
ular we demonstrate that the microcanonical finite-size 
scaling theory presented in Reference for continuous 
phase transitions relates the appearance of these features 
to a peculiar finite-size effect. In addition we deduce 
properties of the dip in the microcanonical entropy which 
we test in numerical simulations of both the Baxter- Wu 
and the four-state Potts model. 

The rest of the article is organized as follows. In Sec- 
tion^a brief introduction to the microcanonical analysis 
of physical systems is given. We particularly discuss the 
properties of the microcanonical specific heat both for 
systems with continuous and discontinuous phase tran- 
sitions in the infinite volume limit. The microcanonical 
finite-size scaling theory of Reference is used in Sec- 
tion mil to deduce the scaling properties of the convex 
intruder in the microcanonical specific heat of systems 
undergoing a continuous phase transition in the thermo- 
dynamic limit. In Section llVl we present numerical results 
both for the Baxter- Wu and for the two-dimensional four- 
state Potts model which belong to the same universality 
class. These numerical results corroborate the findings 
obtained from the finite-size scaling theory as the prop- 
erties of the convex dip are consistent with the theoretical 
predictions. Finally, Section Ivl summarizes our results. 

II. SIGNATURES OF PHASE TRANSITIONS IN 
THE MICROCANONICAL SPECIFIC HEAT 

In finite systems the thermostatic quantities of a sys- 
tem undergoing a phase transition in the thermodynamic 
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limit show characteristic features. The canonical specific 
heat of finite systems, for example, is a regular function 
that exhibits a maximum at a certain temperature. The 
continuous or discontinuous character of the phase tran- 
sition is revealed by the size-dependence of the increase 
of the specific heat maximum which is well understood 
within canonical finite-size scaling theory |Q, . In the 
case of a continuous phase transition the increase of the 
specific heat maximum Cmax with the system size L is 
given by Cmax ^ i"/'' where a describes the divergence 
of the specific heat of the infinite system when approach- 
ing the critical temperature (here and in the following 
we only consider the case a > 0), whereas v is the usual 
correlation length critical exponent. For a discontinuous 
phase transition, however, one finds c,nax ~ L'^ |46l | where 
d is the dimensionality of the system. 

The thermodynamic behaviour can also be investi- 
gated in the microcanonical ensemble which describes 
systems which arc isolated from any environment 0, |3| . 
The starting point in the microcanonical analysis is the 
density of states VI which, for a discrete spin system, mea- 
sures the degeneracy of the macrostate {E), i.e. ^{E) is 
the number of configurations (or microstates) with en- 
ergy E. For a d-dimensional system with linear exten- 
sion L physical quantities are then calculated from the 
microcanonical entropy density 

s{e,L-'):=j-^\nn{L^e,L-^) (1) 

which depends on the energy density e E/L'^ and 
the system size L. In the following the dependence on 
the system size is indicated by the inverse system size 
I :— L~^. Note that the specific entropy also depends 
on the boundary conditions which have to be specified 
for finite systems. This dependence, however, is not in- 
dicated in the notation here. The inverse microcanonical 
temperature shows up as the conjugate variable to the 
natural variable energy e and is defined by 

/3^(e,0 :=^s(e,0. (2) 

The corresponding microcanonical specific heat is given 
by 

c{e,l)^-me,l)r(^^Mllliy\ (3) 

Similar to the canonical case the finite-size behaviour 
of the microcanonical specific heat shows also certain 
characteristics if the system has a phase transition in 
the infinite lattice. In the case of a continuous transition 
the entropy density is concave everywhere. The micro- 
canonical specific heat Q is thus positive and displays 
a maximum whose value increases with growing system 
size L H^ll^. As the derivative of /3^(e, I) is the second 
order derivative of the entropy, it measures the curvature 
of the entropy density. This curvature is negative and has 
a maximum at a certain energy which can be identified as 



the pseudo-critical energy Cpc of the finite system. The 
value of the derivative of (3^ at this pseudo-critical energy 
scales like 

d d^ 

Pfii&pd) ■= -^(3fj,{e,l)\e=epc = ^|j^'S(ej 0|e=epc ^ 

(4) 

in the regime I —>■ 0. The microcanonical critical expo- 
nent Q!^ is thereby related to the canonical critical expo- 
nent a by a, = a/{l -a) El [13. Generally, a canon- 
ical critical exponent k translates to the microcanoni- 
cal exponent = k/(1 — a) if a is positive, whereas 
for non-positive a the critical exponents of the micro- 
canonical description are identical to the canonical ones 
[47l |. Relation is obtained from a decomposition of 
the microcanonical entropy into a regular part and a sin- 
gular part which obeys a certain homogeneity relation 
in the vicinity of the pseudo-critical point. Note that 
the size-dependence of the increase of the specific heat 
maximum allows the determination of the correspond- 
ing microcanonical critical exponent of the specific heat 
of the infinite system. For more details see Section IIIII 
and in particular Reference [2^. Note also that other 
microcanonical quantities like the spontaneous magne- 
tization and the susceptibility also exhibit typical fea- 
tures of continuous phase transitions in finite systems 
El Ei mm m I23. These aspects are however not 
studied in the present work. 

The behaviour of the specific heat is completely differ- 
ent if the system undergoes a discontinuous phase transi- 
tion. In this case the microcanonical entropy density of a 
finite system has a convex dip which leads to a negative 
microcanonic al sp ecific heat for a certain energy interval 
[immilimimill. TWs is different in the canoni- 
cal ensemble where the specific heat, which is related to 
the variance of the energy, is always positive for any finite 
system. Of course in the thermodynamic limit the con- 
vex dip in the microcanonical entropy density disappears 
in systems with finite-range interactions and the specific 
heat is positive everywhere. The curvature properties of 
the entropy density of a finite system allow the deter- 
mination of an inverse pseudo-transition temperature /3t 
by the double-tangent (or Maxwell) construction. This 
is depicted schematically in Figure ^ The discontinuous 
character of the transition is revealed by a certain scaling 
behaviour of the convex dip. Definig the auxiliary func- 
tion s(e,^) := s(e, Z) — /3t (0^7 a size-dependent width cj(Z) 
and depth S{1) of the dip can be measured, as shown in 
Figure^ With increasing system size, i. e. for the regime 
I 0, the depth scales like 6(1) ~ I and the width like 
ti>(/) — ujao ^ I 0, m^ Ell' Here is the length of 
the linear section of the function s of the infinite sys- 
tem and is in fact the latent heat of the transition. The 
asymptotic behaviour linear in I and the appearance of a 
non-zero latent heat ujoo signal the discontinuous charac- 
ter of the transition in the infinite system. This can be 
illustrated by investigating the three-state Potts model 
in three dimensions. The general q-state Potts model is 
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FIG. 1: Schematic depiction of the convex dip in the entropy 
of a finite system whose corresponding infinite system under- 
goes a discontinuous phase transition. The double-tangent 
with slope /3t is indicated by a dashed line (top). For the 
function s — s — Pte the definitions of the width lo and of the 
depth 5 of the dip are shown (bottom). 



defined by the Hamiltonian 



(5) 



where the sum runs over nearest neighbour sites of the 
lattice and the spins ct,; can take on the values ai = 
1, . . . , g |40,|5ll|52l- The coupHng constant is set to unity. 
The three-state model in three dimensions undergoes a 
(relatively weak) first order transition and therefore one 
expects that the characteristic quantities of the convex 
dip evolve linearly in the inverse system size (compare 
|53l | for a canonical and |29l | for a microcanonical inves- 
tigation of the model). This is indeed confirmed by nu- 
merical data as shown in Figure |21 



III. MICROCANONICAL FINITE-SIZE 
SCALING THEORY 

In this Section the implications of the microcanoni- 
cal finite-size scaling theory are investigated for sys- 
tems with a continuous phase transition where > 
(canonically a > v). Note that the condition > 
corresponds to a rather restricted class of model systems. 
Nevertheless important spin models like the Baxter- Wu 
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FIG. 2: The depth Of) and width ujiV) of the convex dip 
in the microcanonical entropy of the three-state Potts model 
in three dimensions as a function of the inverse system size. 
The dashed lines are linear fits to the data. Error bars are 
comparable to the size of the symbols. 



and the four-state Potts model in two dimensions have 
critical exponents cte — and therefore belong to this 
class. These models are closer studied in the next Sec- 
tion. Other models are reported in the literature to have 
exponents satisfying > Ve [s^lssf. 

Consider the entropy density s{e,l) of a finite d- 
dimensional microcanonical system near the pseudo- 
critical energy epc. Recall that epc is defined by the max- 
imum of the second derivative of s. The microcanonical 
finite-size scaling theory for a system with a continuous 
transition assumes that the entropy can be split up 
into two parts, namely into a so-called regular part Sr and 
into a part Sg which is called singular as it becomes non- 
analytic in the infinite volume limit. This is in contrast 
to the regular part which is analytic also in the infinite 
system. The singular part has to satisfy the scaling rela- 
tion 



(6) 



with s :~ e — Cpc and (f> being a scaling function with 
no further explicit size-dependence. This decomposition 
s(£, /) = Si.(e, /) -l- Ss(£, I) is only valid for asymptotically 
large system sizes and in the regime l~^/'^^e — > 0, in 
particular for energies close to the pseudo-critical one. 
Note that the singular part and thus the function (jj in 
relation © is also analytic for finite systems and must 
fulfil the required scaling property which states that it is 
a homogeneous function. Due to the analyticity of the 
microcanonical entropy the decomposition together with 
the required scaling relation for the singular part states 
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that the entropy is of the form 

s{e, I) = so(0 +/3pc(0£ + ImI)£^ + \mI}£^ + ■■■ (7) 

for small e [2^. Here /3pc(0 = Pp-i^pc, 0- The coefficients 
of the second and fourth order terms are of the form 

A2{1) = B2l^ + C2{1) (8) 

and 

Ai{l) ^ B^l^ + Ci{l) (9) 

in terms of the critical exponents of the microcanonical 
system. The coefficient A4 has to be negative, otherwise 
higher order terms in the energy deviation e have to be 
taken into account. The size-dependence contained in the 
coefficients C2 and C4 stems from the regular part Sr of 
the entropy density whereas the dependencies involving 
the critical exponents and have their origin in the 
singular part Sg (or equivalently in the scaling function 
(f) in relation lO) of the entropy. In particular, the coef- 
ficient B2 has to be negative in the case of a continuous 
transition and, due to the regularity of Sr, the coefficient 
C2 is of the form 

C2il) ^ Vil + V2f + . . . (10) 

for small I whereas the coefficient C4 has the expansion 
04(1) = zq + zil + . . . in the inverse system size (see refer- 
ence [2^ for further details). Note that has to satisfy 
< 2 for C4 to be the subdominant term in the coeffi- 
cient A4. This is assumed in the following considerations. 

In certain model systems the critical exponents satisfy 
the relation = Vg. Examples include the Baxter- Wu 
model and the four-state Potts model in two dimensions. 
In this situation the leading behaviour of the coefficient 
A2{1) in relation Q is the linear term in the inverse sys- 
tem size, i. e. ^2(0 ^ ^ foi' ' ~^ where the amplitude a 
is the sum B2+V1. If the amplitude a is positive, i. e. vi 
is positive and large enough compared to the modulus of 
the negative amplitude B2 , the entropy density develops 
a convex dip for all finite system sizes. Note that the cor- 
responding condition A2(/) > is a necessary condition 
for the appearance of a convex dip in the microcanonical 
entropy. Thus, in the framework of the microcanonical 
finite-size scaling theory a convex dip can appear even so 
the transition in the infinite system is continuous. Fur- 
thermore, from expression Q one can then deduce scal- 
ing relations for the width and the depth of the convex 
dip in such a case. Applying the double-tangent con- 
struction to the entropy ((Tjl with = i^^ and a positive 
A2 one finds that 

and 



for small I. This scaling behaviour of the convex dip in 
the microcanonical entropy is therefore different from the 
linear behaviour in the inverse system size encountered 
for discontinuous transitions and discussed in Section UTI 
13 113 • For a negative amplitude a, on the other hand, 
the signatures of a continuous transition with the scal- 
ing behaviour /3^(epc,0 ~ ^2(0 ~ I (= P^/'^O are ob- 
served. Whether the predicted convex dip really appears 
or whether one encounters typical features of continuous 
transitions thus depends on the sign and relative size of 
the amplitudes entering the scaling form {T)) of the en- 
tropy near the pseudo-critical energy. 

Let us end this Section by considering the situation 
where is strictly larger than v^. In this case the 
leading behaviour of the coefficient A2 is still linear in 
I provided the amplitude vi is non-zero. For positive 
vi one then has again a convex dip with the scaling re- 
lations L0{1) ~ Zl/2-(a.-2)/(2^,) g^l^ _ ;2-(a,-2)/>.e 

for its characteristic extensions. For negative vi instead 
one observes features of a continuous transition with the 
asymptotic behaviour /3^ (epc, I) ^ I, which does not allow 
to uncover the actual critical exponents from the leading 
size-dependence of the evolution of the microcanonical 
specific heat, in contrast to the corresponding situation 
for 1^ Ve- 
in conclusion, a convex dip in the microcanonical en- 
tropy is in general possible for systems with a continuous 
phase transition if the critical exponents satisfy > v^. 
For systems with ae > Ve such a dip shows up if the 
amplitude vi of the linear term in the expansion of C2 
in (|10|l is positive. In a situation with ~ the am- 
plitude vi has to satisfy the additional requirement that 
vi > \B2\. In both cases we end up with the necessary 
condition A2{1) > 0, see Eq. ©. 



IV. MICROCANONICAL SPECIFIC HEAT OF 
THE BAXTER- WU AND OF THE FOUR-STATE 
POTTS MODEL 

A. The Baxter- Wu model 

The Baxter- Wu model is a classical spin model whose 
free energy can be calculated exactly in zero magnetic 
field [39! . |40| . The model is defined on a two-dimensional 
triangular lattice, with the Hamiltonian 

n{a) = - ^ a,a,ak (13) 

(ijk) 

where the strength of the coupling constant has been set 
to unity. The spins ai are Ising spins with the possible 
values (Ji = ±1. The sum in H13|) extends over all trian- 
gles of the triangular lattice (here denoted by (ijk) with 
i,j and k specifying the vertices of the triangles). The 
system has a continuous phase transition in the infinite 
system at the critical temperature Tc = 2/ln(l -I- ^/2). 
The critical exponents which characterize the singular 
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behaviour of physical quantities near the critical point 
are known exactly. For the specific heat and the cor- 
relation length one has for example a = v = 2/3 or 
= = 2 microcanonically. In this Section the micro- 
canonical specific heat of the Baxter- Wu model is inves- 
tigated numerically [H^. To this end the microcanonical 
entropies of systems with linear sizes ranging from 16 
to 96 are calculated using the very efficient transition 
observable method [23, llj for both open and periodic 
boundary conditions. 

As the Baxter- Wu model exhibits a continuous phase 
transition in the infinite lattice, one expects at first sight 
that the microcanonical specific heat shows a maximum 
with the associated increase for increasing system sizes. 
For finite systems with o'pen boundary conditions this is 
indeed the case. The derivative /3^(epc,/) of the micro- 
canonical temperature at its maximum value, i. e. at the 
pseudo-critical energy, is negative and its modulus de- 
creases, leading to an increasing specific heat maximum 
when the system size is increased. For the Baxter- Wu 
model the critical exponents (Xe and are equal and 
therefore one expects from (0J and the discussion in Sec- 
tion mil to observe the scaling relation /5^(epc,0 ~ ^ fo'^ 
small I. This expectation is indeed confirmed by the data 
shown in Figure |31 
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FIG. 3: The derivative /3^epc,0 = d/3^(epc, 0/de of the 
Bcixter-Wu model as a function of the inverse system size I for 
open boundary conditions. The errors are approximately as 
large as the symbol size. The dashed line shows a linear fit. 
The data extrapolate to zero in the limit I 0, corresponding 
to a divergent specific heat in the infinite system. 



The situation is different for finite Baxter- Wu systems 
with 'periodic boundary conditions. In Figure^the func- 
tion s is shown for a finite system with linear extension 
L = 60. A convex dip is clearly visible. This observation 
apparently suggests at first sight that the model seems 
to undergo a discontinuous transition in contradiction to 
the known continuous character of the phase transition. 
Note that recently a double-peak structure in canonical 
histograms of the energy was reported in the literature 

El. 

In the previous Section IIIII we have shown that the 
microcanonical finite-size scaling theory permits the ap- 
pearance of such a convex dip for the case = I'e- With 
the relations H12|) and (|ll|l and the known value = 2 




FIG. 4: The function s for the finite Baxter- Wu system with 
L = 60 and periodic boundary conditions (here Pt = 0.43945). 
Note that s is only obtained from a numerical simulation up to 
an additive constant. The data show a double-peak structure 
typical of finite systems with a discontinuous transition in the 
macroscopic limit. 



for the Baxter- Wu model we expect to observe the scal- 
ing behaviour ~ for the width and 8{l) ~ P for 
the depth of the convex dip. These predictions are in- 
deed confirmed by our numerical data obtained for the 
Baxter- Wu model as shown in FigureEl The full lines are 
fitted functions containing the expected asymptotic be- 
haviour in the inverse system size I. For the dashed lines 
the leading correction term coming from the regular part 
Sr of the microcanonical entropy (which is proportional 
to /"^ for the depth and to p/'^ for the width) has been 
included. In agreement with the continuous character of 
the phase transition the dip scales away in the macro- 
scopic limit I Q. The size-dependences of both the 
depth and the width of the convex dip follow the theo- 
retical predictions. This is especially clear from the inset 
of Figure [S] where we compare lo{1) with the expected 
behaviour in a double- logarithmic plot. This agreement 
gives strong evidence to the general scaling theory pre- 
sented in the previous Section ITTTI 

It should be stressed again that the different behaviour 
of the microcanonical entropy of the Baxter- Wu model 
for different boundary conditions is related to different 
properties of the expansion coefficients in Q). In partic- 
ular the coefficient A2 seems to be deeply affected by the 
boundary conditions. A microscopic explanation of this 
observation is however still lacking. 

The findings of this section have direct implications for 
the behaviour of canonical energy histograms, which have 
been discussed for the Baxter- Wu model in Reference 
m. The energy histograms are directly related to the 
canonical distribution function Ppie) which is the proba- 
bility of finding a microstate with energy density e at the 
inverse canonical temperature p. In terms of the micro- 
canonical entropy this distribution function is basically 
given by Pp{e) ^ exp (/3L''e — L'^s{e, I)) . Thus for the in- 
verse pseudo-transition temperature Pt obtained from the 
Maxwell construction applied to the microcanonical en- 
tropy the canonical energy histogram exhibits two peaks 
at equal height with the peaks appearing at the same 
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FIG. 5: The depth S{1) and width uj{l) as functions of the 
inverse system size I for Baxter-Wu systems with periodic 
boundary conditions. The errors are of the order of the sym- 
bol size. The full lines show fits to the theoretically predicted 
behaviour (IIH and 11211 using the six largest system sizes. 
The dashed lines include in addition the first correction term 
coming from the regular part Sr of the microcanonical en- 
tropy. The inset shows a double-logarithmic plot of the width 
together with a straight line of slope 1/2. From the six largest 
systems one gets a slope of 0.481(3). 



energies as the peaks in s = s — /?te. It follows that 
in the canonical histograms the peaks at the energies ei 
and 62 approach each other according to the scaling law 
|ei - 62! ~ VI. 

An important quantity which influences the kinetics 
of first order phase transitions is the interface tension E 
defined by 



2/3tS = 



1 



■In 



^'/3t(emax) 
^'/3t(emin) 



(14) 



where the double-peaked distribution Pp^ has one of 
the two peaks of equal height at energy Cmax whereas 
the minimum between the two peaks appears at the 
energy Cmin [13 • From the above discussion one gets 
S ~ ~ / and thus the interface tension vanishes in 

the limit of asymptotically large systems. Consequently 
a coexistence between an ordered phase and a disordered 
phase at temperature /3t is hardly detectable in this case 
although the canonical distribution function Pp^ exhibits 
a double-peak structure. For moderately small systems, 
however, the dip represents still a considerable barrier be- 
tween the order phase and the disordered one. In some 
sense one can speak of a phase coexistence for not too 
large finite systems. 



B. The four-state Potts model 

The four-state Potts model in two dimensions belongs 
to the same universality class as the Baxter-Wu model. 
The character of the phase transition of the (7-state Potts 
model in two dimensions changes exactly at g = 4 from 
continuous for low values of q to discontinuous for large 
q. This makes the four-state Potts model particularly in- 
teresting. Due to the presence of marginal scaling fields 
for g = 4 the singular behaviour of the four-state Potts 
model is not simply described by pure power laws, but 
physical quantities like the specific heat or the suscepti- 
bility acquire additional logarithmic factors [58ll59l |. Log- 
arithmic corrections then also show up in the finite-size 
scaling form of the specific heat • This is in contrast 
to the Baxter-Wu model where these marginality effects 
are absent and the singular behaviour of the specific heat 
does not exhibit logarithmic corrections [53, |60| • 
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FIG. 6: The derivative /3^j(epc, /) = d/9^(epc, i)/de of the four- 
state Potts model as a function of the inverse system size I 
for open boundary conditions. The errors are approximately 
of the order of the symbol size. The dashed line shows a fit 
with an asymptotic linear term and corrections including a 
logarithmic term. The data extrapolate to zero for the limit 
/ — > 0, yielding a divergent specific heat in the infinite system. 
The inset shows a fit without the logarithmic correction term 
(see main text for further details). 



Using the same numerical technique as for the Baxter- 
Wu model, we have computed the microcanonical en- 
tropy density of the four-state Potts model on a square 
lattice with open and periodic boundary conditions. For 
open boundary conditions system sizes ranging from 
L = 12 to L = 80 have been considered, whereas for 
periodic boundary conditions system sizes from L = 14 
to L = 50 have been simulated. 

As shown in Figure El the derivative /3^(epc,0 of the 
microcanonical temperature evaluated at the pseudo- 
critical energy is found to be negative for open boundary 
conditions. As its modulus decreases as a function of the 
inverse system size, the resulting specific heat maximum 
increases when the system size is increased, thus yielding 
the typical finite-size features of a continuous transition 
in the thermodynamic limit. The critical exponents 
and i^e being equal for the four-state Potts model, one ex- 
pects (see Eq. (0J) the scaling relation /?' (epc,/) ~ I for 
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small I. However, as Figure demonstrates, this asymp- 
totic regime is not yet fully reached for the system sizes 
considered here. Presumably, this has its origin in the 
presence of logarithmic correction terms. Nevertheless, 
the data are consistent with the expectation of a lead- 
ing linear dependence in the inverse system size in the 
asymptotic limit / 0. To demonstrate this a function 
of the form (/3|i(epc, 0)fit — al + bl/\n{l) + cP has been 
fitted to the data (dashed line in Figure 0). To get a 
rough impression of the presence of a logarithmic correc- 
tion, a fit function with correction terms biP + cil^, i. e. 
without a logarithmic term but with the same number 
of fit parameters, has also been fitted to the numerical 
data, yielding a much poorer result (sec inset of Figure 

El). 

For the Potts model with periodic boundary conditions 
our findings are in qualitative agreement with the con- 
siderations of Section IIIII and with the results obtained 
for the Baxter- Wu model with periodic boundaries. The 
microcanonical entropy displays the same convex dip for 
finite systems. Note that a double-peak structure in en- 
ergy histograms obtained from canonical Monte Carlo 
simulations was already reported in the literature for this 
model ■ The scaling behaviours of the depth and the 
width of this dip are again consistent with the predicted 
scaling relations and H12() as depicted in Figure 
even so we were forced to include a logarithmic correction 
term in addition to the leading correction term coming 
from the regular part. 
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FIG. 7: The depth S{1) and width uj{l) as functions of the 
inverse system size I for four-state Potts systems with periodic 
boundary conditions. The errors are approximately as large 
as the symbol size. The dashed lines show fits of the expected 
scaling relations lllll and il'Zt to the data with a logarithmic 
correction term. 



V. CONCLUSIONS 

The determination of the order of a phase transition 
presents a basic issue in the statistical analysis of phys- 
ical systems. As there are only few model systems that 
can be tackled analytically this question is usually ad- 
dressed through numerical simulations. A widely used 
indicator of a first order transition is the metastability 
between two coexisting phases [s^ [s^, ^3- This leads 
to a convex dip in the microcanonical entropy and thus 
to a negative microcanonical specific heat in the transi- 
tion region. In the canonical ensemble the convex dip 
corresponds to a double-peak structure in histograms 
of suitable observables. However, the appearance of a 
double-peak structure in canonical histograms or a con- 
vex dip in the microcanonical entropy is not sufficient 
to identify a first order transition as such features have 
been reported in the literature for certain model systems 
that are known to undergo a continuous transition in the 
thermodynamic limit. To identify the first order charac- 
ter additional scaling properties of the convex dip in the 
microcanonical entro py ( or the double-peak structure in 
canonical histograms jS^ll^l) have to be satisfied. 

In this work we analyzed the Baxter- Wu model and the 
two-dimensional four-state Potts model, which are known 
to have continuous transitions in the thermodynamic 
limit, but for which double- pea k structures in canoni- 
cal histograms are observed |4ll . li^ . Using a recently 
developed phenomenological finite-size scaling theory for 
systems with continuous phase transitions |2y| the pos- 
sible properties of the microcanonical entropy of mod- 
els where a > v have been investigated. It turns out 
that finite-size effects may indeed lead to a convex dip 
in the entropy of any finite system although the mod- 
els show continuous transitions. The finite-size scaling 
relations for the microcanonical entropy predict in this 
case certain scaling properties of the convex dip that are 
different from those encountered at a first order transi- 
tion. The numerical studies of the Baxter- Wu system and 
the four-state Potts model with periodic boundary con- 
ditions indeed reveal an entropy with a convex dip. The 
properties of the convex dip are in agreement with the 
expected scaling relations deduced from the microcanon- 
ical finite-size scaling theory. This gives strong s upp ort 
to the validity of the scaling theory proposed in [2y| for 
finite microcanonical systems. The systems with open 
boundaries on the other hand do not show a convex dip 
in the corresponding entropy. For this case too, the com- 
puted specific heat satisfies the scaling relation of the 
microcanonical scaling theory. The different behaviour 
for different boundary conditions could be traced back 
to the qualitatively different properties of the expansion 
coefficient A2 in (7J). Future studies have to clarify in 
more details the influence of different boundary condi- 
tions on the finite-size behaviour of the microcanonical 
entropy function. 

Even so the appearance of a convex dip in the micro- 
canonical entropy usually signals a first order transition. 
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it may also show up in systems with a continuous transi- 
tion. We have shown in this work that its appearance in 
those cases (and therefore also the related appearance of 
a double-peak structure in canonical histograms) can be 
traced back to a finite-size effect and that the properties 



of the dip, which are markedly different from those as- 
sociated with a first order transition, can be understood 
within a scaling theory developed for the description of 
continuous phase transitions in the microcanonical en- 
semble. 
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